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A semi- classical description of the intrinsic spin-Hall effect (SHE) is presented which is relevant for 
a wide class of systems. A heuristic model for the SHE is developed, starting with a fully quantum 
mechanical treatment, from which we construct an intuitive expression for the spin-Hall current and 
conductivity. Our method makes transparent the physical mechanism which drives the effect, and 
unifies the SHE across several spintronic and optical systems. Finally, we propose an analogous 
effect in bilayer graphene. 
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Introduction. The spin-Hall effects^ (SHE) entail a col- 
lection of phenomena in which a pure, transverse spin 
current is generated in response to a longitudinally ap- 
plied electric field. Unlike extrinsic versions of the effect^ 
the intrinsic SHEs are finite in the absence of impurities, 
originating from the spin-orbit coupling (SOC) present in 
the band structure of the host system. The rich study of 
the SHE has been largely fueled by its potential applica- 
tion as a spin current source in the emergent technology 
called spintronics.^ 

The intrinsic SHE has been studied in many systems, 
especially in semiconductor spintronic systems. Mu- 
rakami et al. studied the SHE of holes in p-doped bulk 
semiconductors.^ Independently, Sinova et al. studied 
the effect in two-dimensional electron gases (2DEG) with 
Rashba SOC.^ Remarkably, the spin-Hall conductivity 
(SHC) in this system was found to be universal, i.e. 
independent of system parameters, which subsequently 
stimulated much research of the SHE in 2DEGs (having, 
in general, both Rashba and linear Dresselhaus SOC)^^!^ 
The cubic Dresselhaus SOC case was studied by Bernevig 
and Zhang. ^ Two-dimensional hole systems with Rashba 
SOC also exhibit the intrinsic SHE,^^^^'^^ a unique signa- 
ture being a resonant SHC.^^ The SHC is typically char- 
acterized by the linear Kubo response of the transverse 
spin current to a charge current excitation. However, this 
treatment conceals the physical mechanism which drives 
the effect. The SHE in bulk ^-semiconductors^ for exam- 
ple, arises from the Berry curvature of momentum space 
which results in spin-dependent equations of motion. 
Similar in nature are the recently proposed SHEs of 
light, phonons,^^ and excitons.^^ On the other hand, 
the SHE in Rashba 2DEGs occurs as a result of time- 
resolved spin dynamics induced by an electric field. We 
focus in this paper on the latter type, and hereafter we 
shall use "SHE" to mean this type. 

In this paper, we start with a general spin-orbit model 
and provide a semi-classical picture of the SHE, explain- 
ing clearly how it is driven. From there, we construct 
an intuitive expression for the spin-Hall current and con- 



ductivity. Our general expression reproduces exactly the 
SHC obtained previously using linear response theory in 
a wide class of systems. Hence our approach not only 
provides a clear, physical picture for the SHE, but also 
allows one to classify various SHEs as originating from a 
common underlying mechanism. Furthermore, our closed 
expression for the SHC can conveniently be applied to 
other systems using simple vector algebra, without the 
need for cumbersome operator algebra. 

Theory. We consider the general spin-orbit Hamilto- 
nian in the presence of an electric field. 
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where p = hk is the momentum, m the effective mass, 
7 is the SOC strength, a is the vector of Pauli matri- 
ces, Ct{k) is a momentum-dependent effective field, and 
E is the electric field. We study the time {t) evolution 
of the above quantum system. To incorporate the ex- 
plicit t-dependence of the system quantum mechanically, 
we switch to the interaction picture, splitting H into 
two parts, H = Ho -\-Hi, where Ho = eE - f governs the 
time evolution of the operators, and Hi = |^ —ja-Cl{k) 
governs the time evolution of the states in the new pic- 
ture. The momentum operator in the interaction picture 
(subscript /) is found to be pi{t) = e^'^ot/h^^-iHot/h _ 

p — eEt^ i.e. with the expected linear t-dependence due 
to E. State vectors |'0(t)) in the Schrodinger picture 
correspondingly transform as |'0/(t)) = e*^°^/^|?/^(t)), 
and evolve according to the new "Schrodinger equation" , 
Hi (t) lip I {t)) = ihdt\il)i{t))^ where the Hamiltonian Hi{t) 
is found to be 



Hiit) = ^-ja-n(ki(t)), where 
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and the summation over repeated indices is implied. The 



2 



Hamiltonian Hi{t) is that of a particle subject to an ex- 
plicitly t-dependent spin-orbit field, Q{t). We consider 
diagonalizing the Schrodinger equation at time t, by ap- 
plying a unitary rotation U{t) which rotates the z-axis 
to be pointing along the instantaneous field i.e. 

U{t)Hi{t)U\t) = U{t){ihdt)U\t), 
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where e = ihdt is the energy operator. On the right- 
hand-side, a gauge field Ao{t) = —iUU^ appears from 
the t-dependence of U. The term hAo{t) represents an 
additional Zeeman-like term, indicating the presence of 
an effective magnetic field in the rotating frame. Since 
U = ^iUa ' uj [18], where uj is the instantaneous angular 
velocity of the rotating frame, the effective magnetic field 
in the laboratory frame is found to be —Co. Clearly, the 
angular velocity of the rotating frame is also that of the 
unit vector n{t) = Ct{t)/\Q{t)\^ so = n x n. Thus, the 
effective magnetic field arising from the gauge transfor- 
mation is Cl± = nx which is an invariant with respect 
to the particular choice of /7, depending only on the t- 
dependence of Q{t). We illustrate this field in Fig. [H 
A classical derivation of this field can be found in Refs. 
mil^. The Q± term represents a physical magnetic field 
which couples to the electron spins, ^^^^^ and, as we show 
below, is precisely the component leading to the SHE. 
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FIG. 1: (Color online) In the presence of a time-dependent 
magnetic field, B{t) = |5(t)|n(t), an additional magnetic field 
B± = ft X n (green, vertical arrow) is seen by spins. The net 
instantaneous magnetic field felt by spins is the vector sum of 
B{t) and B±, denoted by the dashed, black arrow. 

Semi- classical calculation of spin- Hall current. We 
consider the DC response of a SOC system to a charge 
current ji due to an electric field Ei. In particular, we 
calculate the spin current j^j = s/2{vj^a^} where s is the 
value of the spin angular momentum, and also the SHC 
alj = jlj/Ei. The electric field results in a t-dependent 

spin-orbit field Ct{t)^ which, as we showed above, is al- 
ways accompanied by an additional field Q± = ft x ft 
(Fig. [1]). We assume that the spins align themselves to 
the net magnetic field (Ij], i.e. the sum of Q and the t- 
dependent correction. The classical spin vector is then 
s = :^sCIj:/\CIj:\^ where +(— ) represents spins aligned 
parallel (anti-parallel) to Ctj^. Along the /-coordinate, 
the component of the spin vector due to ^}_\_ ii 



±s (^n X where is expressed in terms of its 

equivalent angular velocity. We now consider the adia- 
batic limit. In the ideal limit, oo, so the spins 

always remain aligned to it as it varies. In reality, \Q\ is 
finite and there is a secondary component and the 
relevant condition is ^ i.e. the spins are pri- 

marily aligned to yet with a small component along 
\Ct±\. The adiabatic spin polarization along / is 

I 

^ , ^, o ^n^r^lnr (3) 

where einr is the Levi-Civita symbol. We employ the con- 
ventional spin current operator, j^ j = s/2{vj^ a^}, where 
{•} denotes the anticommutator. The velocity operator 
reads Vj = dH/dpj = pj/m — jh~^a • dCl/dkj^ and upon 
substituting this expression into the spin current opera- 
tor we obtain j^ j = s (pja^ /m — jdQi/dkj) . The spin 
current therefore has two contributions, the first due to 
the spin-polarization of the states, and second due to the 
variation in k of Cl{k). The spin-Hall (sH) current is the 
part of jl j proportional to Ei^ which corresponds only 
to the first contribution, thus ,^ (sH) = pjS^/m. The to- 
tal spin-Hall current can be calculated semi-classically by 
summing the expectation value of the j^ ^ (sH) operator 
over all states up to the Fermi level. 



A,w = /^£f),'(*,, (4) 

where D is the dimension of the system. We carefully 
note that the expression for the spin polarization s^{k) 
has opposite signs for spins pointing parallel and antipar- 
allel to Ct{k) ([3|). These correspond to the two spin- 
split bands of the spin-orbit Hamiltonian ([Tj); namely, 
electrons pointing parallel (antiparallel) to Q have en- 
ergies ^+(-) = p^ /2m — (+)7|1](^)|. Since the Fermi 
level is usually exceeds the band-splitting, Ep ^ 2\Q.\^ 
both bands are occupied with their Fermi wavevectors 
fulfilling Ef = {hk$)'^ /2m±j\n{k$)\, where k^ > k%. 
Thus, in the region where the two Fermi surfaces overlap 
(i.e. k < /c^), there is a complete cancellation of s^{k). 
The finite contribution to the spin-Hall current comes 
from the states in the Fermi surface occupied only by 
the ground state band, E_. Noting this, and writing the 
time derivative of in the expression for (|3|) in terms 
of the wavevector h, we obtain 
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Dividing both sides through by the electric field Ei, we 



where the integration in k goes from limits k^ n.^. 



obtain an expression for the intrinsic SHC, 
(sH) esh f kj 
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This is our central result. Using a semi-classical ap- 
proach, we have constructed a general expression for 
the intrinsic SHC which makes transparent the physical 
mechanism driving the phenomena. Our analysis is rel- 
evant for describing the SHE in a wide class of systems, 
as shown below. 

Results and Discussions. In Tab. we list several sys- 
tems which can be represented by the general Hamilto- 
nian in Eq. ([1]). In each case, we give the expression 
for the ^-dependent magnetic field Q{k)^ and compute 
the i-spin polarization, s^{k), using Eq. (|3]) assuming an 
electric field applied along x. Since the s^{k) are odd in 
ky^ carriers traveling in opposite ^-directions become po- 
larized along opposite i-directions. This is the physical 
picture of the SHE. Below we determine the SHC in each 
of the spintronic systems. We also examine the optical 
analog, and finally propose an analogous effect in bilayer 
graphene. 

A. Combined Rashba and Dress elhaus SOC — The out- 
of-plane spin polarization driven by an electric field in the 
combined Rashba-Dresselhaus system is given in Tab. . 
It vanishes when \a\ = in agreement with previous 
studies.^ When |a| 7^ using Eq. (|6|), and the relation 
k- = -2m\nRD\/h^, the SHC is 
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in agreement with Ref. [7]. In particular, for the a 
or /3 = cases, the universal value is reproduced. 

B. n-doped bulk semiconductors — For conduction 
electrons under the influence of fc^-Dresselhaus SOC, the 
SHC © is 
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where is the /c^-Dresselhaus field, and (O^cj)) are 
spherical angles in ^-space. Using the interband relation 
k^ — k^ ^ —2m\^{kF)\/h'^kF, where kp = {k^ -h kp)/2, 
we produce the SHC of cr^y = kF/127r'^ as computed in 
Ref. 0. 

C. Holes in UI-V semiconductor quantum wells with 
Rashba SOC — The SHC ([6l) is (5 = 3/2), 
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where n is the hole density, giving a universal SHC of 
^xy ~ —9e/S7r}^ In comparison to the value of —e/Sir 
in 2DEGs, the extra factor of 9 here arises from the 
larger angular momentum of heavy holes, and the k^- 
dependence of the spin-orbit coupling. 

D. Rayleigh scattering of polaritons — An analogous 
effect occurs in optics, when polaritons undergo Rayleigh 
scattering. The polariton polarization is represented by 
a pseudospin r , where the pseudospin field is exactly that 
of bilayer graphene (see below). Upon scattering (which 
changes the wavevector /c), polaritons acquire a finite 
component, corresponding to circular polarization, whose 
sign depends on the initial k. This is the optical SHE*^ 

E. Bilayer graphene — Finally, we propose an analo- 
gous effect in bilayer graphene (BG). The BG system is 
modeled as two coupled honeycomb lattices, with each 
layer having two inequivalent lattice sites A, B and A, B 
in the top and bottom layers respectively. We assume 
the Bernal stacking {A-B) configuration. In the low en- 
ergy limit, electrons in the BG system are described by 
an effective 2-by-2 Hamiltonian,^^ Hbg = ' ^bg-, 
where bg = {k'^ — ky^2kxky^0). The r here is the vector 
of Pauli operators acting on the pseudospin, rather than 
the actual electron spin. The BG eigenstates ^ are two 
component wavefunctions describing the electronic am- 
plitude of electrons on the two layers, ^ = ('^(A), tlj{B)). 
When an electric field E = E^x is applied to the bilayer 
system, an additional out-of-plane component accompa- 
nies the strictly in-plane pseudospin field CIbg- The 
out-of-plane component induces in the adiabatic limit a 
pseudospin polarization along z (see Tab. ), which corre- 
sponds to charge transfer between the two layers. This 
is an essential ingredient for the technology known as 
pseudospintronics^^^ in which binary states are encoded 
by relative charge densities on the two BG layers. The 




FIG. 2: (Color online) Illustration of proposed pseudospin- 
Hall effect in bilayer graphene. The arrows indicate the di- 
rection of the electron momenta, (left) With no electric field, 
electrons with all momenta are distributed evenly between 
the two layers, (right) With an applied electric field in the 
x-direction, electrons are separated to each layer depending 
on their ^-momentum; electrons with -\-{—)py > are trans- 
(9) ferred to the bottom (top) layers respectively. 



where (j) is the azimuthal angle in ^-space. In the limit of 
small Rashba coupling, the band resolved Fermi wavevec- 
tors can be shown to be k^ ^ ^J2i\n ^ (2Am/fi^)7rn, 



physical effect of this polarization is illustrated in Fig. [2j 
electrons with > (j)y < 0) are separated to the bot- 
tom (top) graphene monolayer of the BG system. The 
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TABLE I: List of systems and their Hamiltonians Ti in which the spin-Hall effect is analyzed. D is the system dimension, 
a\l — z) are the Pauli spin operators, ki are the wavevectors, (t± — ^ ia'^ and k± = kx ^ iky. Q{k) is the momentum- 
dependent effective magnetic field, and is the z-spin polarization of electrons resulting from an electric field applied in the 
x-direction, obtained using Eq. (|3|). For the case of bilayer graphene, we mean here the pseudospin polarization, . 



effect is completely analogous to the SHE with replace- 
ments a ^ r. It may therefore be called the pseudospin- 
Hall effect Such an effect should be of interest to the 
field of pseudospintronics. 

We finally reiterate that there exists a second class of 



the intrinsic SHE arising from the non-trivial curvature of 
^-space, which was not discussed in this paper. ^^ ^^^i^^i^^ 
We have shown elsewhere that the underlying mecha- 
nisms, although distinct, are related. 
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Although distinct, the two mechanism are actually related 

The spin may well have a component along / due to Q itself, 
but it is not proportional to the electric field and does not 



contribute to the spin-Hall current. 



